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Abstract 

We study here what it means for the Universe to be nearly flat, as 
opposed to exactly flat. We give three deflnitions of nearly flat, based on 
density, geometry and dynamics; all three definitions are equivalent and 
depend on a single constant flatness parameter e that quantifles the notion 
of nearly flat. Observations can only place an upper limit on e, and always 
allow the possibility that the Universe is infinite with = — 1 or finite with 
k = 1. We use current observational data to obtain a numerical upper 
limit on the flatness parameter and discuss its implications, in particular 
the "naturalness" of the nearly flat Universe. 

Keywords: general relativity — cosmology 

1 Introduction 

Observational cosmologists tell us that, at present, the ratio of the sum of the 
densities of all forms of matter-energy in the Universe to the critical density 
required for spatial flatness is: 

f^T.o = 1.02 ±0.04 . (1) 

We follow standard usage throughout, using the symbol H. to denote densities 
relative to the critical density p<-rit = (3/87rG')(d/a)^, with Pa — A/SttG repre- 
senting the contribution from vacuum energy so that 17^ = (p + Pk)/ Pciit- The 
subscript "0" indicates terms taken at the present time. 

Eq. 1^ is a 95% confidence-level fit to combined data on fluctuations in the 
spectrum of cosmic microwave background (CMB) radiation from the WMAP 
satellite and observations of the magnitudes of distant Type la supernovae 
(SNIa), plus the assumption that the present value of Hubble's parameter satis- 
fies Ho = {a/a)g > 50 km s^^ Mpc^"'^ 1 . The data are consistent with the pos- 
sibility that r^T.o is precisely one, implying an infinite and spatially flat {k = 0) 
Universe. They are also marginally consistent with the idea, widely held during 
the 1990s, that Sl^.o < Ij implying a Universe which is infinite and negatively 
curved (fc = —1). Most intriguingly, Eq. ^ suggests the possibility of a closed, 
finite and positively curved Universe with fJ^.o > 1 and k = +1 — a prospect 
which has received less theoretical attention over the years but which 
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Figure 1: Closed (fc — +1), flat (fc = 0) and open (k — —1) universes. 



is now hinted at by both the location of the primary acoustic peak |H1 Ej and 
the amplitude of low-order multipoles (HI El in the CMB power spectrum. These 
are three radically different models of the Universe (Fig. 1). 

The flat fc = Universe is widely favored among cosmologists for two rea- 
sons. First, inflation dynamically causes a very small portion of the primordial 
Universe to grow so large as to encompass the entire observable Universe, so that 
the spatial structure of the Universe is greatly flattened. This, however, only 
forces the Universe to be nearly, rather than exactly flat ^-^d indeed 

there are examples of inflationary models with fc 7^ ^1 El El El 
Secondly, it is felt by some that a nearly flat Universe involves a numerical co- 
incidence or fine-tuning between 17^ and unity, especially at early times. As we 
will discuss in §6, however, this apparent flne-tuning is a natural result of the 
definitions of the density ratio JIt and the critical density, combined with the 
cosmological equations. We must also mention that notions of naturalness and 
fine-tuning are subjective; for example, one might consider a flat Universe to be 
infinitely fine-tuned since it has Six identically equal to one, thereby making it 
the most unnatural choice. 

Our main purpose in this work is to define just what it means for the Universe 
to be nearly fiat. We give three logically independent but physically equivalent 
definitions, based on density, geometry and dynamical behavior. All lead to the 
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same dimensionless constant, whose significance as a flatness parameter does 
not appear to have been widely appreciated. ExpHcitly, this parameter is the 
ratio of the de Sitter radius of the Universe to a constant of integration in the 
Friedmann-Lemaitre equation. 

We derive our flatness parameter on the basis of density in §3, geometry in 
§4, and dynamics in §5. The question of naturalness is addressed in §6, and in 
§7 we give our summary and conclusions. 

Our discussion has features in common with that of Chernin Chernin 
discusses four parameters, which he calls Friedmann constants; these are Ay 
associated with vacuum energy, Ajj associated with dark matter, Ab associated 
with baryons, and Aji associated with radiation. He emphasizes that all are very 
roughly equal (within a few orders of magnitude), and suggests that this may 
be evidence of a cosmic symmetry of unknown origin. On the other hand, we 
emphasize that for a nearly flat Universe, Ay must be at least several orders of 
magnitude less than Ad, and that indeed this is responsible for the near flatness. 
Thus our work is consistent with and complementary to that of Chernin. 

The issue of near-flatness is also raised in a recent article by Lake who 
adopts a rather different approach from ours, based on the trajectories of model 
universes in the fi^ — fi^ plane. Lake's intent is to assess the likelihood that 
the Universe is exactly flat, whereas we would contend that observation cannot 
distinguish — even in principle — between a perfectly flat Universe and one 
that is sufficiently close to flat. 



2 Cosmological Dynamics 

We begin by briefly reviewing cosmological dynamics, focusing on the present 
epoch in which the Universe is dominated by a cosmological constant (or dark 
energy) and pressureless matter. The matter component could be further broken 
down into cold dark matter and baryonic matter, but these have the same effect 
on the dynamics and we treat them together. For early times, radiation is also 
important jl8j . We assume a linear equation of state for the matter component, 
p = wp with w = const. The basic equations of homogeneous and isotropic 
cosmology may be written (with c = 1): 

SttGp = -a + 3 f 4 + 4 V 8^Gp = A-(JL + ^ + ^] , (2) 

\a a J \a'^ a a J 

With the linear equation of state and some algebraic manipulation, we obtain 
a first integral of these (the Friedmann-Lemaitre equation) as 
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Figure 2: Scale factor as a function of time in the standard concordance model. 



where C is a constant and is the de Sitter radius, defined by: 

3 ^ _J_ ^ _J_ ^ 1 , . 

A SttGpa h^^a H^nA.o ■ ^ ' 

Here r^A.o — Pa/ Pciit.o is the present density ratio of vacuum energy; we note 
that the vacuum energy density is constant but the density ratio JIa is not. 
(The constant C here is the same as in Lake's paper 19 , while C and Ra are 
referred to as and Ay respectively in the paper by Chernin ^Hl-) Most of the 
discussion in this paper will be based on Eq. ©. For the case of a flat (fc = 0) 
Universe, it can be solved exactly when the matter component is dominated by, 
e.g., pressureless dust {w — 0) or radiation (u> = 1/3): 

- / sinh2/=^(3t/2i?A) {w = 0) 

(Ci?A)^/^sinhi/2(2i/i?A) (^«-l/3) ■ 

The first of these expressions governs the now-standard concordance model of 
present-era cosmology (see Fig. 2). It is found in surprisingly few of the standard 
texts EOlEIl and goes back to Heckmann in 1931 [^1^ . 

Some special times are of particular interest in this model: the time when 
acceleration begins, the time of equality between energy densities of matter and 
dark energy, and the present time (when dark-energy density is approximately 
2.8 times the matter density). These are conveniently characterized by: 

{I/V2 (acceleration begins) 

1 (equal matter/dark-energy density) . (6) 

yJ^A,o/^m,a (at present, « 1.67) 

Best- fit values of the relevant cosmological parameters from the WMAP satellite 
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Figure 3: Behavior of the function /(a) that determines the density ratio, with 
minimum = i{C/2Rj,f/^ at a = [cRl/2f/^. We have used Eqs. Q and 
1)12(1 for i?A and C respectively, with WMAP values for f^M.o and JIm. fl- 
are Ho — 71 km s^^ Mpc^^, r^M.o — 0.27 and r^A.o — 0.75 P, so that the de Sitter 
radius works out to i?A — 16 Gly from Eq. Q). 

Note that for the case fc = the equations are scale-invariant, so the scale 
function is arbitrary to within a multiplicative factor and is not a measurable 
quantity. For k ~ ±1 the scale function is measurable in terms of H and 0,^, as 
we will show in the next section. 



3 Density Definition of Nearly Flat 

Since density can never be measured with perfect precision, it is clear that 
we cannot use the quantity Qt.o itself to verify whether or not the Universe 
is precisely flat. But the observational bound Q certainly indicates that the 
Universe is either flat, or nearly so. In fact, we will use fi^ to obtain our first 
definition of nearly flat. 

The first of the cosmological field equations (0) together with the first inte- 
gral (jSJ allows us to solve for 0,^. in terms of the scale factor as follows: 

r k 

(C/a)3-+i+a2/i?2 

where the function /(a) reads 

f{a) ^ (C/a)3-+i ^ ^2/^2 ^ ^-/a + / RI (8) 

for the present era with w — Q. We plot this function in Fig. 3 and note that its 
minimum value is 3(C/2i?A)2/'^. The total density parameter itself is plotted in 



1 - 



/(«) 



(7) 
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Figure 4: Behavior of the density ratio versus scale factor, using Eqs. Q and 
H12() for i?A and C respectively, and assuming WMAP density parameters at 
present [a/ao = 1): JIt o = 1-02 ± 0.04. The flat case JIt.o = 1 is also shown for 
comparison. 



Fig. 4; note that 17^ — > 1 in the early Universe when a — > 0, and that 17^ ^ 1 in 
the distant future when a becomes very large. Since observations indicate that 
r^T is close to one, the second term in square brackets must be small, so we may 
expand as 

The parameter Vt^. will be close to unity if 

/(a) = C/a + a^/Rl > = 3(C/2i?^)2/3 » i . (10) 

A nearly flat Universe with a total density parameter Q^- that is nearly equal to 
unity for all times is therefore equivalent to one with a small flatness parameter 
£, defined as 

£EEi?A/C<l. (11) 

This parameter is the constant ratio of a fundamental constant in the theory 
(Ra) to a constant of integration (C), so it is allowed by the theory to have any 
value. (It is equal to Chernin's Av/Ajj 18 .) However, many theorists would 
consider it unnatural for e to differ from unity by many orders of magnitude, or 
for C to differ by many orders of magnitude from the natural value Ra- 

The constant C is simply related to the matter density; this can be seen by 
comparing the first of Eqs. (|2Jl with Eq. ^ divided by (for w = 0): 

C = {8TTG/3)ipa^) = n^H\^ = n^.oH^al, (12) 
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where the final steps follow from the definition of and the fact that C = con- 
stant. This relation expresses the well-known way in which matter density de- 
creases as the Universe expands. 

Let us assume that the Universe is not exactly flat, but is nearly flat, and 
use the observational bound to obtain an upper limit and rough estimate 
for the flatness parameter. If we divide the first of Eqs. by 3ff^ we obtain 

- 1 = klH^a^ . (13) 

Combining our expressions I0J and p2|) with this, wc find that 

e = [Hlal^^^.^f^X^ = < 0.012 . (14) 

It then follows from the definitions ^ and (|ll|l that the constant C must be at 
least 

C = ^= > 1300 Gly . (15) 

Moreover the scale function can be determined explicitly if fc 7^ 0; from Eq. H13|l 
it is given by 

k 1 



which takes the value Oq « 99 Gly at present. (This result nicely dislays how 
a becomes indeterminate for fc = and fl-r = 1.) The value of the flatness pa- 
rameter in Eq. H14|l is small, but can hardly be considered unnatural. Similarly 
comments apply to the distances in Eq. ((T^ and (fTH|l . 

It is amusing to note that an independent means of determining e is available 
in principle. From the first of Eqs. we may obtain 

iH/H„f = r!M,o(ao/a)' + ^a.o - {^t.o - l)ia„/af . (17) 

For k = 1, {H/HoY dips below its asymptotic value of f^A oi dropping to a 
minimum of {H.JH„f = [1 - [4,/21)e^]n^^„ at {a^/ao) = in^,j2{n^^„ - 1). 
Thus observations of the expansion rate at this epoch would yield a definitive 
measurement of e. Unfortunately the relevant time occurs so far in the future 
(of order 60 Gyr) that such a method is not of practical interest. 

In summary, we may say that, in terms of density, the Universe is nearly 
flat if the flatness parameter deflned in Eq. (|11|1 is small, which means that 
the constant of integration C is large in comparison to the de Sitter radius R^- 
Future improvements in measurements of the density ratio J^t.o are obviously of 
great interest, particularly since its current value hints at a density higher than 
the critical one, implying that the Universe may have positive spatial curvature 
and be finite. 
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4 Geometrical Definition of Nearly Flat 



Geometry provides the most intuitive method for defining near-flatness. A part 
of the surface of a sphere is nearly flat if its extent is smaU compared to the 
radius of the sphere. In the same way, we may say that the Universe is nearly 
flat if the region to which we have access is small compared to its characteristic 
radius. For a metric in the standard form, 

ds^ = de - a" (t) ( + u^de^ + sm^e d(p'\ , (18) 
\ 1 — ku'^ J 

this may be expressed mathematically by saying that the Universe is flat within 
some region (i.e., within proper radius (.) if its dimensionless comoving radial 
coordinate u remains small throughout that region: 



u 



< 1 . (19) 



Our region of interest is the accessible Universe, so we take t to be the Hubble 
distance I = c/H. 

Proper (or physical) distance ^ is calculated in terms of u by 



a{t) du' 



a{t) sin "^u (k = 1) 
a{t)u (k = 0) 



V 1 - ku'^ [ a(ijsinh"u (fc = -1) 

^ a{t)u[l + |u2 + . . . ] . (20) 



Eq. (|20|l is easily inverted for u = u{£): 



u 



= {e/a)[l-l{e/ar + ■■■]. (21) 



We will need only the lowest order in {£/a), so that the value of k is irrelevant. 
Taking £ = 1/H = a/d for the accessible Universe, we flnd for the square of the 
comoving Hubble distance: 

ul = (1/df - {C/a + a^/Rl - k)-^ = [/(a) - k]-^ . (22) 

The equivalence of the definitions in terms of density and geometry becomes 
clear if we use Eqs. JTl and to relate to follows: 

< - ^ . (23) 

That is, the comoving Hubble distance and the deviation of the density ratio 
from unity are small together, and both vanish at very early and very late times. 
Using the WMAP value fi^.o ~ 1.02, the present spatial extent of the Universe 
in terms of the comoving Hubble distance is Uh,o ~ 701)2 = 0.14. 

In brief summary, the geometric definition of near-flatness is completely 
equivalent to the density definition. 
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Figure 5: The mechanical analogy for the dynamics of the Universe for each of 
the three cases of spatial curvature. The maximum of the potential occurs at 
X = 1/21/3 where V;„„ = -3/22/3. 

5 Dynamical Definition of Nearly Flat 

The first integral ^ may be put into dimensionless form by the substitutions 
T = t/Rj, and x = a/{CRiy/^ to obtain 

{dx/drf -{l/x + x^) = -fce2/3 . (24) 

It is apparent that, for small e, the dynamical behavior will be similar to that 
for fc = 0. Thus, as in the previous sections, we are led to e as a measure of 
near-flatness. Eq. can of course be solved exactly by quadratures, as we will 
do below, but for a perturbative solution the appropriate expansion parameter 
is clearly ke^/^ . 

It is worth noting the mechanical analogy for Eq. (|24() . which is a particle 
of mass 2 and total energy E = —ke^^^ in a potential V{x) — —{1/ x + x'^), as 
shown in Fig. The particle flies rapidly outward from the origin, coasts over 
the top of the potential hill, and accelerates down the hill to infinity. Obviously, 
if the total energy E is much smaller than the distance between the top of the 
hill to zero (which is 3/2'^/3), then total energy is unimportant and the three 
cases /c = 0, ±1 will exhibit the same qualitative behavior. The difference is a 
shift in the time required to reach a given value of x. 

Let us illustrate the use of Eq. 124|l . Very early in the history of the Universe, 
its expansion is dominated by the \/x term, while for its very late history it is 
dominated by the x'^ term. Only at intermediate times is the last perturbative 
term relevant. This is easily seen if we take the positive root of Eq. I|24|l and 



9 



solve exactly by quadratures, then expand in S = fce^/^ as follows: 

Jo y^l/x + X^-. 

The integrals /„ are given by 



J2 ^"Cn/„ . (25) 



^" - / (1/y + j'2)„+l/2 ' (26) 



with cq — 1 for n ~ and c„ = ^ 2-4-'^(2n) '' ^'^^ n > 0. 
To zeroth order in S one has 



r = /o = (2/3) ln(.T3/2 ^ + ^3) . (27) 

This is equivalent to the first of Eqs. lO, valid for fc = 0. The higher-order 
terms in the expansion (|25|l represent the shift in time required for the Universe 
to attain a given value of a when k = ±1, relative to the fc = case. The 
integrals J„ cannot be expressed in terms of elementary functions for n > 0. 
However, a change of variables to s = allows one to rewrite them in terms of 
the hypergeometric function F{a, b; c, z), following Silbergleit {24) : 

/„ = i / (l + s)-("+l/2)s"/3-l/2^^ 



i?g^F(n+i,t + i;t + |,-x3). (28) 



In particular, it is possible to solve exactly for the accumulated time shift very 
late in the history of the Universe using Eq. H26(l . In the limit as a goes to 
infinity, one has: 

hm / - r - + (29) 

— ""io (i/y + yT+1/2 - 3r(n + i) ■ ^^^^ 

The leading contributions to At = t — Iq are then found (in the late-time limit) 
as 

lim At = i(fce2/3) lim + |(fce2/3)2 + ■ ■ ■ 

= 0.29 fce^/^ + 0.075 k'^s*^^ ■ ■ ■ (30) 

In conjunction with the upper limit (|14|l on e, the contributions of these two 
terms to the asymptotic shift in physical time At = RaAt amount to no more 
than ±240 Myr and 3 Myr respectively. 

The time shift At that has accumulated up to the present time can be 
calculated by evaluating Eq. H26(l numerically with x — Xg — {fl/^ o/^MfiY^^', one 
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finds in this way that the first two terms in the expansion contribute no more 
than ±120 Myr and 2 Myr, respectively. These times are small on cosmological 
scales, but not negligible. For instance, they are much longer than the duration 
of the radiation-dominated era. 

6 Flatness of the Early Universe 

At present the Universe is either flat with fc = 0, or nearly flat with flatness 
parameter e ^ 10~^. But the early Universe was extremely close to density- 
flat, as emphasized by many authors, notably Dicke [25) . Some authors have 
suggested that fc = is the most natural choice, since otherwise the density 
ratio JIt must have been fine-tuned to be very close to unity, to within about 
a part in 10^"^ at the time of nucleosynthesis. This is the "flatness problem" 
of standard cosmology, and one of the issues that the inflationary paradigm is 
claimed to address [TU] . 

Here we advocate a rather different view, for two reasons. First, as mentioned 
in the introduction, the notion of naturalness is fundamentally subjective, and 
one is equally entitled to view the case fc = as infinitely fine-tuned, with 

identically equal to one. Second, as we will show below, the difference 
between the actual density and the critical density actually diverges like in 
the early radiation-dominated Universe. The density ratio fix approaches unity 
only because the critical density diverges even faster, like t~^. The apparent fine- 
tuning is thus inherent in the definition of the critical density and the equations 
of general relativity, in particular Eq. Q. As we will discuss in detail below, 
flrj. must approach very close to one at early times. Thus the quantity fi^ — 1 is 
not a useful measure of naturalness or fine-tuning, and we suggest the use of the 
constant fiatness parameter e instead. Since this is as large as 10~^ according 
to present observations, a nearly fiat Universe cannot as yet be claimed to be 
excessively fine-tuned. 

Let us first calculate the deviation of the total density from critical in the 
early Universe. From the first of the cosmological equations (0) and the defini- 
tion of the critical density, we see that 

= (3/8^G)(fc/a2) . (31) 

In the early radiation era, a cx t^^^ so for k ^ the difference in Eq. H31|l 
diverges at early times like t~^. The quantity — 1 = {pt — P^iit)/ P„it is small 
only because the critical density diverges even more rapidly: 

= (3/87rG)(a/a)2 cx _ (32) 

The smallness of the quantity f^T — 1 at early times is thus an artifact of the 
definition of the critical density. 
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Let us next estimate the deviation of the density ratio from unity at the 
time of decouphng, which occurred at a redshift of about « [ao/a^) ~ 10^ at 
a time t « 10^^ s; this is also roughly the time of matter and radiation density 
equality. From Eq. Q: 

- 1 « k/f{a) « ka/C . (33) 

This is a reasonable approximation in the matter-dominated era (i.e., until re- 
cently), during which a oc t^/'^ . Thus the deviation of from unity at decou- 
pling obeys 

(17^,, - l)/(r!^,o - 1) « flaK « « 10^3 . (34) 
Since the present deviation is about ^ 10~^ we obtain a small deviation of about 
< 10~^ at decoupling. 

Similarly, we may go back in time in the radiation era, during which w = 1/3 
and a cx t^/'^, to a time of about t„ « 1 s, when nucleosynthesis occurred. As 
above, but with w = 1/3 in Eq. iQ, we obtain 

(r!T.„ - l)/(f^T.d -l)~ {aja^f « tJU ~ 10"'" . (35) 

Thus the deviation during nucleosynthesis was only about 10~^'', an impressively 
small number. It follows merely from the present density deviation, which is 
not necessarily very small, plus the cosmological equations and the definition of 
the critical density. 

A useful analogy from elementary physics might be the following: consider a 
test particle of mass m with total energy E falling into the Newtonian gravita- 
tional field of a mass M . The ratio of this particle's kinetic energy K = mv^/2 
to its potential energy |[/| = GMm/r is K/\U\ = {E/GMm)r + 1. Note that 
the difference K/\U\ — 1 becomes arbitrarily small as one approaches r ^ 0, 
in exactly the same way that fi^ — 1 does in cosmology as t — > 0. Yet one 
would hardly be justified in concluding from this that E "must be" zero on the 
grounds of naturalness. 

In summary, the extremely small deviation of the density ratio from unity in 
the early Universe is a consequence of the definition of the critical density and 
the basic equations of relativistic cosmology for any value of k. We therefore 
do not agree with the viewpoint that fc = is necessarily the most natural 
interpretation of current observational data. If future experiments produce a 
much smaller limit on the flatness parameter e (say, 10~^), then that might be 
a more convincing indication that the most natural value for k is zero. 

Chernin has made similar observations, with which we concur although 
our interpretation is somewhat different. Chernin emphasizes that the most 
notable feature of the Universe is not its near- flatness, nor the densities of its 
various components, but the closeness of the Friedmann constants C and i?A to 
each other. He has subsequently argued that this apparent coincidence may in 
turn be understood most naturally in a closed universe whose radius of curvature 
is of the same order of magnitude as its present size |2fij . 



12 



7 Summary and Conclusion 

We have presented a new quantitative definition of a nearly-flat Universe in 
terms of a flatness parameter e, and shown that it follows uniquely from three 
independent lines of argument based on density, geometry and dynamical be- 
havior. It is clear from our derivations that there is no way, even in principle, to 
distinguish between a precisely flat Universe and one whose flatness parameter 
is sufficiently small. Measurements of total density from CMB data currently 
set an upper limit on e of order 10^^, which is small but not extremely small 
or unnatural. By contrast, the variable quantity VL^- — 1 is not a useful flatness 
criterion, because JIt is necessarily driven toward unity at both early and late 
times (for any value of k) by the equations of relativistic cosmology and the 
deflnition of the critical density. 
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